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Abstract 



A mathematical model to describe period-memorizing behavior in Physarum Plas- 
modium are reported. In constructing the model, we first examine the basic charac- 
teristics required for the class of models, then create a minimal linear model to fulfill 
these requirements. We also propose two modifications of the minimal model, non- 
linearization and noise addition, which improve the reproducibility of experimental 
evidences. Differences in the mechanisms and in the reproducibility of experiments 
between our models and the previous models are discussed. 
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1 Introduction 



The plasmodium of Physarum polycephalum is a huge unicellular organism 
with amoeboid movement. It is differentiated into an advancing front zone and 
a rear region that is composed of a network of protoplasmic veins. T he concen- 



trations of many chemical com ponen ts in its cytoplasm, such a s Ca 2 + (IRidgwav et al 



19761 : lYoshimoto et allll981ah. ATP flYoshimoto et allll981bL H+ flNakamura et al 



1982h . cAMP dUeda et allll986h . and NADH flMori et allll987h . display mutu- 



ally entrained oscillations in a period of 1-2 minutes. These oscill ations also in- 

volve the contraction-relaxation cycle of cytoplasmic actomyosin tlWohlfahrt-Bottermannl . 
19791 ). which generates rhythmic shrinkage of protoplasmic veins and the ac- 
tive transport of cytoplasm via shuttle streaming in the vein network. The 
streaming of cytoplasm pushes the front zone to spread out, and the solation 
of cytoplasmic gel in the zone occurs at the same time; thus the plasmodium 
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migrates. This shuttle transport of cytoplasm also coordinates the phase re- 
lations among local oscillatory d ynamics and lets a hug e plasmodium behave 
as a unified individual organism (I Yoshimoto et all Il978l ). 



Because of its particular behaviors, the plasmodium of P. polycephalum has 
been intensively examined as a model organism for the study of the coop- 
erative phenomena in biological systems. Many theoretical approaches have 
also been reported, i n which the p l asmodium is modeled as sp atially extended 



oscillatory med ium (ITeplov et all Il99ll; lYamada et al. 



dilat or system flAkahane et all Il999l : iTero et all l2005t lOdelll . 11982 ; ISmith 



1999) or coupled os- 



19941 ) . A combination of the results of experimental and theoretical stud- 



ies have revealed a lot of interesting features abou t Plasmodium ' s beha viors 



such as entrain ment to environmental s timulation tjMiyake et all Il992l). mu 



tual entrainment flYoshimoto et all 1 1978T). symmetry-induced patterns (ITakamatsu et al 



2001 
2008h 



Takamatsul . 120061 ) . maze-solving (lNakagaki et all 120001 ). and computing tlAono et al. 



Recently, another interesting behavio r - memorizat i on of t he period of periodic 
stimulation - has been reported by Saigusa et al.l ( 2008 ). In that study, the 
Physarum plasmodia were exposed to dry and cold condition (dry stimulation) 
three times periodically with periods T = 30-90 minutes. The stimulation re- 
duced the locomotion speeds of the Physarum plasmodia. After the last stimu- 
lations, spontaneous slowdowns of locomotion speeds were observed. The time 
intervals between the last stimulations and spontaneous slowdowns coincided 
with the stimulation periods. This phenomenon was referred to as a spon- 
taneous in-phase slowdown (SPS). The number of SPS occurrences variated 
from none to three, depending on the organism. Additionally at T = 60min 
after SPS stopped, one-time stimulation induced the spontaneous slowdowns 
after an equal period of time from the additional stimulation, which is referred 
to as SPS after one disappearance (SPSD). SPSD was not observed at T = 40 
min or T = 80 min. These experiments clearly indicate that the Physarum 
plasmodium memorizes the period of the periodic stimulation. Furthermore, 
since the dry cold conditions used in the study are thought to be unfavorable 
for Physarum plasmodium, the spontaneous slowdow n can be interpreted as 
anticipated behaviors against unfavorable conditions (ISaigusa et all 120081 ) . 



In this paper, we present mathematical models to explain this period-memorizing 
behavior. Since there is little information about the controlling mechanisms 
of the Physarum Plasmodium's locomotion, we assume a general dynamical 
system, and examine the basic characteristics required to control the Plasmod- 
ium's locomotion speed. We then present a simple model that fulfills these 
requirements and qualitatively reproduces experimental evidences. We insist 
the model is minimal in the sense that it has the minimum degrees of free- 
dom and that the functions that it uses (including the differential equations 
to drive the system) are all linear functions. We also propose some modifi- 
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cations of the minimal model to improve the reproducibility of experimental 
evidence. Another mat hematical model has been presented in a paper report- 



ing these experiments (jSaigusa et all 120081 ) . though that model's mechanism 



differs from ours. We outline the previous model and compare it with our 
model in the Discussion section. 



2 Basic requirements for models 



In this section, we argue the class of models for a biochemical system which 
controls the reduction of the locomotion speed of the Physarum Plasmodium. 
The main characteristics of the system are that the stimulation makes the 
system reduce the locomotion speed, and the system also spontaneously re- 
duces the speed. Given that the speed reduction is a transient response, the 
systems behaviors are analogous to those of a neuron with steady and puls- 
ing inpu ts. Thus , our m odel refers to simple neuron models: integrate-and-fire 



models (IBurkittl . 120061 ) . We consider differential equations 



±i = g(x u x 2 , ■ ■ ■ ,x n ), i = l,---,n, (1) 

to describe the autonomous dynamics of the biochemical system in P. poly- 
cephalum, and suppose that this system reduces the locomotion speed when 
the variables satisfy a condition 



f{x 1 ,x 2 , ■ ■ ■ ,x n ) > 0. (2) 

In other words, f(x%, x 2 , ■ ■ ■ , x n ) = determines the threshold plane, and if the 
system gets across the plane, locomotion speed is reduced. If the stimulation is 
imposed, the orbit of the system in the sub-threshold region is forced to shift 
into the over-threshold region. On the other hand, if the orbit spontaneously 
crosses over the the threshold plane, SPS occurs. 

In the following, we examine the experimental evidences, list four basic re- 
quirements for models, and give appropriate assumptions to meet these re- 
quirements. 

Equilibrium states. Before the periodic stimulation was applied or long after 
stimulation was ceased, the Plasmodium showed steady locomotion. Thus, the 
system should have equilibrium states in the sub-threshold region (f(xi, x 2 , • • • , x n ) < 
0), so that it does not reduce speed spontaneously in the equilibrium states. 

Locomotion recovery. Physarum plasmodium recovered its locomotion speed 
after stimulation, and became ready to respond to the next stimulation. This 
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implies that the condition eq. ([2]) does not hold for long after the slowdown 
begins. In order to bring the system's state back to the sub-threshold re- 
gion, here we assume that some of the variables are reset after the sys- 
tem induces the slowdown. Let x\, • • • ,x m be the variables to be reset and 
Hi = Xm+i, ■ ■ ■ , Vn-m = %n the variables not to be reset after the slowdown. 

Transitive spontaneous action. SPS did not persist in the experiments. 
Thus, there should be a difference in the system's actions between a stimulation- 
induced slowdown and spontaneous slowdown (SPS). We assume that the re- 
setting manners differ. In the SPS case, the reset may start before the orbit 
invade deeply into the over-threshold region. However, the stimulation can 
be vigorous and may shift variables widely so that the function ([2]) takes a 
large positive value. This value may affect the resetting manner. In particular, 
we assume that after the stimulation, resettable variables are shifted to the 
complete reset state 



reset _ . 1 /0 n 

Xi — > x u i = 1, • • • ,m. (3) 

On the other hand, the reset after SPS is incomplete: resettable variables are 
shifted to an internally dividing state between the state just before the reset 
and the complete reset state 



Xi r ^> XiXi + (1 - Xi)Xi, i = 1, • • • , m. (4) 



Memory. The Plasmodium remembered the stimulation period for a long time 
after stimulation stopped. In order for a dynamical system to memorize such an 
analog value, it must have a slow manifold in the phase space: the relaxation 
along with the manifold is quite slow so that the syst e m can preserve the 



information about its own history for a long time (iFoxl . let all ). As a simple 
and extreme case, we introduce the neutral manifold (a line of equilibrium 
states) in the phase space of our model, where the system never relaxes in 
that direction, and thus, the stored memory is permanently preserved. In 
other words, we introduce a conserved quantity whose value never changes in 
the autonomous dynamics. Along with those dynamics, the system's orbit is 
constrained on a n — 1 dimensional invariant plane transverse to the neutral 
manifold, and the resetting action will shift the orbit to another plane along 
with the neutral manifold. The period of stimulations is thought to control 
which plane is selected. 

Summarizing the class of models, we can say that these models have a conser- 
vative quantity; thus the autonomous dynamics of the models are restricted 
to the n — 1 dimensional invariant plane. In this plane there is an equilib- 
rium state, a n — 2 dimensional threshold plane, and a complete reset state. 
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Changing the value of the conservative quantity changes the positional rela- 
tion among them in the plane. The relation determines whether or not au- 
tonomous dynamics from the complete reset state to the equilibrium state 
cross the threshold plane, and if they do, it also determines the period of time 
from the complete reset state to the threshold plane. 



3 Models and their behaviors 



Based on the above discussions, here we present a minimal model for the 
slowdown induction system. The model has two resettable variables x±, X2, 
and one un-resettable variables y, driven by linear differential equations, 



Xi — x 2 - Xi 

t x x 2 = y - x 2 (5) 
Ty y = Xl -y 

where r x and r y indicate the relative time scales of x 2 , and y respectively. 
These equations have a conserved quantity 

X\ + t x x 2 + r y y = constant. (6) 

Thus there are neutral manifold in the phase space, and the equilibrium states 
form a line 

x 1= x 2 = y. (7) 
We set the threshold function as 

f(x 1 ,x 2 ,y)=x 1 -y-5, (8) 

where S determines the distance between the equilibrium state and the thresh- 
old plane. We choose the complete reset state for resettable variables 

(x 1 ,x 2 ) = (0,l). (9) 
Thus, reset after stimulation becomes 



Xl ^0 

(10) 

reset _, 
X 2 > 1 
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and reset after SPS becomes 




reset' 



reset' 



Ai ■ xi 



A 2 • x 2 + (1 — A 2 ) 



(11) 



We show the behavior of our model in fig. [TJ Initially the system stays at 
the equilibrium point x\ = x 2 = y = 1. We give the stimulation three times. 
After that, at t = 4.03, the system reaches the point where eq. (jSJ) changes 
the sign from negative to positive, SPS occurs, and the orbit relaxes to a new 
equilibrium point x\ = x 2 = y — 0.51. After additional stimulation is imposed 
(t = 10), SPSD occurs at t — 11.26. and the orbit relaxes to an equilibrium 
point that is closer to the previous one (x\ = x 2 = y — 0.52). 

We repeated the simulations with various stimulation periods. The response 
periods of SPS and SPSD are shown in fig. [2J The parameters used in these 
simulations are tuned so that these response periods have similar values to the 
stimulation periods and show positive correlations as the stimulation period 
changes. 

Here we sketch the spontaneous dynamics of the system, and verify that the 
above model with three variables is a minimal model. Because the system 
has one conserved quantity, the spontaneous dynamics are restricted in a two 
dimensional plane, which is described in fig. [3J The plane is divided by thresh- 
old line eq. (jSJ) and has one equilibrium point and one complete reset state 
which lies in the sub-threshold region. Initially the arrangement of the plane 
may look like fig. [3]A.. The relaxation orbit from the complete reset state di- 
rectly relaxes to the equilibrium point. Periodic stimulation may change the 
positional relationship among the equilibrium point, the complete reset state, 
and the sub-threshold line, similar to fig. O-B. That occurs due to the con- 
tinuous change of the unresettable variable during the periodic stimulation. 
In the plane, the orbit from the complete reset state violates the threshold 
line in the course of relaxation to the equilibrium point, and thus SPS occurs. 
After the induction, the orbit shifts a bit toward the the complete reset state 
and relaxes to the equilibrium point. We emphasize that after periodic stim- 
ulation, the orbit from the complete reset state overreaches and crosses the 
threshold line, but the orbit after spontaneous induction directly relaxes to 
the equilibrium point. This happens because the plane is two-dimensional. If it 
is one-dimensional, there are only two cases: if the equilibrium point is in the 
sub-threshold region, the orbit relaxing to it never crosses the threshold, and 
if the equilibrium point is in the over-threshold region (this case is eliminated 
in the above assumptions), the orbit never reaches the equilibrium point. This 
ensures that the plane should be at least two-dimensional. Thus, with the 
neutral manifold for memory, a minimal model must have three variables. 
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To make the model minimal and easily understandable, we chose linear func- 
tions for differential equation eq. dSJ), the threshold function eq. (|Hj) and re- 
setting manner eq. (II ip . However, these functions can be nonlinear without 
changing the basic mechanisms. Here we give an example of the nonlinear 
modification of the model, by letting parameters r x , r y depend on variables 



11. . 
► — • (1 + oj x1 ■ x x + u x2 ■ x 2 + uj x3 ■ y) 

T f T f (12) 

► — ■ (1 + U y i • X X + U y2 • X 2 + UJy3 ■ y) 

and changing the threshold function from eq. (jHJ) to 



fniin{xi, x 2 , y) = xi - y ■ (1 + ugi ■ %i + uj S2 ■ x 2 + uj sz ■ y) - 8 (13) 

where Wi are the nonlinearity parameters. The modification improves the re- 
producibility of the periods, as displayed in fig. [21 

In experiments, the number of times of SPS occurrences varied among; whereas 
one organism showed SPS three times, another organism did not at all. These 
variations must be due to fluctuation in the cell state dynamics. Adding 
Langevin noise to the differential equations 



Xi = x 2 - Xi + cr£i(£) 
r x x 2 = y - x 2 + a£ 3 (t) (14) 
r y y = x x -j/ + erf 3 (t) 

where (^(^1)^(^2)) = % ' ^(^1 — h), our model also shows such varieties of 
responses, as shown in fig. HI 



4 discussion 



In this paper, we presented a mathematical model for the period-memorizing 
behavior of the plasmodium of P. polycephalum. We described and exam- 
ined the conditions for satisfying the model, and introduced a minimal linear 
model that qualitatively reproduces the experimental evidence. To improve 
the model, we also suggests the two methods (nonlinearization and noise ad- 
dition) that are better for quantitative reproducibility. 

Since almost nothing has been revealed about the chemical components for the 
locomotion speed reduction system of the Physarum plasmodium, we did not 
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present a detailed model in which variables correspond to actual chemical com- 
ponents. Instead, we made a minimal linear model whose mechanism is easily 
comprehensible. In our model, for example, the system's resetting action after 
it reduces the locomotion speed is modeled in a conditional branching man- 
ner, instead of in an autonomous dynamical system. In a detailed model, this 
process may be represented by neuron-like firing action. The abstraction revel 
of our model c orresponds to those of integrate-and-fire models for a neuron 
( IBurkittl . 120061 ). as ment ioned abo v e. The transcriptome of P. polycephalum 



is now being sequenced (iGlochnerl . let all ). This information may lead us to 
develop a more detailed model. In particular, two important actions for lo- 
comotion, actomyosin contraction and sol-gel transition of the cytoplasm, are 
thought to be controlled by the calcium ion concentration. Thus, transcripts 
homologous to what are involved in calcium ion concentration regulations in 
other organisms may give useful insight into the detailed mechanism. 



Previously, anot her model for the Pla smodium's period-memorizing behavior 
was presented by lSaigusa et al.l (120081 ) . Here we outline that model's architec- 
ture and behaviors. That previous model consists of decoupled multiple oscilla- 
tors. The oscillators form sub-groups, each of which contains a certain number 
of oscillators with the same frequency (oc 1/period). There are multiple sub- 
groups with different frequencies. The phase relations among oscillators are 
random at the initial condition. The periodic stimulation first synchronizes 
oscillators in phase within sub-groups that have similar periods to those in 
the stimulation, and then it synchronizes these sub-groups in phase. In that 
study, this super cluster (a cluster of clustered sub-groups) is assumed to in- 
duce SPS in their study. After stimulation is ceased, the super cluster gradually 
desynchronizes and SPSs are observed as only transient phenomena. However, 
the clusters within the sub-groups do not desynchronize, because they have 
exactly the same periods. An additional one-time stimulation again induces 
the formation of a super cluster, and thus induce SPSD. In the following, we 
remark on the two advantages of our model in reproduction of experimen- 
tal evidence, and discuss two other prominent differences in the mechanism 
between the previous model and our model. 



With the previous model, the occurrence of SPSD depends on at which phase 
the one-time stimulation is given; stimulation at the in-phase with previous 
periodic stimulation induces SPSD much more strongly than the stimulation 
at the anti-phase. However, this is not consistent with the experimental results. 
In the experiments, the occurrence of SPSD did not depend on a small shift 
in the timing of one-time stimulation 0, which is reproduced in our model. 
Besides, SPSD was observed only at T = 60 min in the experiments, and 
our model also gives SPSD in a certain window of the stimulation period 



the occurrence of SPSD seems to decrease monotonically when the delay in timing 
increases. This may be explained by the noise effect in both models. 
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(T = 0.84 ~ 1.06 with the linear model and T = 0.82 ~ 0.98 with the 
nonlinear model). However, the occurrence of SPSD in the previous model 
does not depend on the stimulation period. Thus another mechanism will be 
required for the previous model to reproduce the occurrence of SPSD for a 
certain window of the stimulation period. 

There are large differences in required degrees of freedom, between our model 
and the previous model. As described above, a response period in the previous 
model is represented by intrinsic collective modes of the dynamical system. 
This means that a number of variables (oscillators) are needed to respond to 
each stimulation period. Thus the model needs a lot of degrees of freedom. 
Actually, 4.4 x 10 5 os cillators are used in the numerical simulation in the paper 



( jSaigusa et all 120081 ): it contains 440 sub-groups with different frequencies and 
each sub-group consists of 1000 oscillators. In our model, on the other hand, 
a response period is not generated by the intrinsic mode but is represented 
by the length in phase space between the reset state and the threshold plane, 
which is continuously adjustable. Thus our model shows qualitatively the same 
behaviors with only three degrees of freedom. 

Because the previous model uses intrinsic modes, the stimulation period and 
response period show fairly good correspondence. No tuning of parameters is 
necessary for the correspondence. On the other hand, in our model, parameter 
tuning is necessary even for reproducing the qualitative behaviors. It should 
be noted that having a conserved quantity in a general dynamical system itself 
demands the tuning o f parameters, since the system is structurally unstable 



( 1 Guckenheimerl . let al.l ). However, if this period-memorizing qualification has 
an adaptive role as means for anticipating unfavorable conditions, it might be 
not so hard for biological organisms to tune these parameters in evolution. For 
example, the period and the response to the temperature changes of circadian 
oscillation sys t ems a re thought to be tuned evolutionarily in many organisms 



f lRosato et all 119971 ). 
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Fig. 1. Time series of the minimal linear model. 

(T x ,T y ,S,X 1 ,X 2 ) = (1.11,1.23,0.0961,0.571,0.592) are used. Red line indi- 
cates the time series of x±, green line indicates X2, and blue line indicates y' = y + 5 
instead of y (so that spontaneous crossing over of the threshold is identified by 
crossing of red and blue lines). Stimulations are induced at t = 1,2, 3, 10. Arrows 
with dashed lines indicate the reset of resettable variables after the stimulations, 
and arrows with solid lines indicate the reset by the spontaneous crossing over the 
threshold. 
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Fig. 2. Response periods after various periods of stimulations. Data from 
minimal linear model and from nonlinear model (in which eq. (|12p and (|13p 
are used instead of eq. ([5]) and ([XT]) ) are plotted with red and blue re- 
spectively, with parameters (t x , T y , 8, Ai, A2) = (1.11,1.23,0.0961,0.571,0.592) 
for the linear model and {t' x , r y , 6', Ai, A 2) u x x, uj x2 , u X 3^yl, ^y2, ^y3^Sl, ^S2, ^63) = 
(0.926, 1.32, 0.0672, 0.46, 0.805, 0.116, -0.111, -0.0803, -0.118, 0.151, -0.089, 0.047, 0.0603, 0.0577) 
for the nonlinear model. The lengths of time between the last stimulation of peri- 
odic stimulation and SPS are plotted with filled squares, and the lengths of time 
between the one time stimulation and and SPSD are plotted with filled circles. 




Fig. 3. The schematic figures of the invariant planes transverse to the neutral 
manifold and orbits on them. Plane A (displayed at left) corresponds to initial 
unstimulated conditions, and plane B (right) corresponds to the conditions after the 
periodic stimulation ceases. The white area indicates the sub-threshold region, and 
the shaded area indicates the over-threshold region. The black open circle indicates 
the complete reset state, and the black filled circle indicates the equilibrium state. 
The red arrow shows the relaxation orbit from the complete reset state to the 
equilibrium state. In the plane B, the orbit crosses the threshold line (the subsequent 
trajectory is drawn by dashed line). The orbit then shifts a bit toward the complete 
reset state, and relaxes to the equilibrium state. Both are drawn by blue arrows. 
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Fig. 4. Number of SPS responses in the linear model with Langevin noise, eq. 
(PT3l). (r x ,Ty,S,Xi,X 2 ) = (1.11,1.23,0.1361,0.571,0.592) and a = 0.04 are used for 
the simulation. For each stimulation period, 100 simulation runs are performed and 
statistics on the number of SPS responses are taken. The magenta area represents 
the proportion of cases showing three SPS responses, blue shows the proportion of 
cases with two SPS responses, and green and white show the single-response and 
no-response proportions, respectively. 
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